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Considering i/e — s- and i/e ^ i^fj. oscillations as the signals in neutrino factory we study the 
discovery reach of CP violation in presence of standard and non-standard interactions of neutrinos 
with matter. For standrad interactions it is found that around 3000 Km baseline the discovery 
reach of CP violation is better for neutrino factory. Even this is better than the 130 Km baseline 
in superbeam facility but with the conservative choices of neutrino flux and detector specifications 
as given by GLoBES. However, the precision measurement of CP violating phase S in neutrino 
factory is found to be good in the baselines ranging from 3000 Km to 4000 Km. In presence of 
non-standard interactions the discovery reach of CP violation for 3000 Km baseline is better than 
130 Km superbeam facility in presence of noon-standard interactions and also for most of the 
allowed region of Est- Ifowever, 130 Km baseline is found to be better for e^^ and e^^- For other 
NSIs - See , £tt the 3000 Km baseline in neutrino factory is better for smaller values of NSIs whereas 
130 Km baseline with superbeam facility is better for larger values of NSIs. Compared to other 
baselines for neutrino factory the NSI discovery reach is in general better for 3000 Km baseline 
except for e^r where as for example, 2300 Km baseline for superbeam is found to be better. 
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I. INTRODUCTION 



The present experiments on neutrino oscillations confirms that there is mixing between different flavours of 
neutrinos (vei v^)- The probability of neutrino oscillations depends on various parameters of the neutrino 
mixing matrix-the PMNS matrix [1]. The current experiments tells us about two of the angles 623 and 612 
[2] with some accuracy. The reactor neutrino experiments like Daya Bay[3] and Reno[4] provided compelling 
evidences for a relatively large angle ^13, with 5.2a and 4.9(7 results respectively. These recent reactor 
neutrino results indicate 613 very close to 8.8°. The CP violating phase 6 is totally unknown. Although the 
mass squared difference of the different neutrinos {Amfj — mf — m'j) are known to us but the sign of Am§i 
(which is related to mass hierarchy) is still unknown. 

In this work we consider long baselines {L > O (2000) Km) for neutrinos coming from neutrino factory. 

However, it is difficult to consider short baseline below 2000 Km for nciitrino factory as the first oscillation 
peak of oscillation probability will correspond to energy of order MeV for which proper neutrino flux will not 
be available. We consider totally active scintillator detector (TASD) for our study. Considering — > i^n and 
i^e I'm oscillations as the signals we do a comparative study of discovery reach of CP violation for different 
long baselines in neutrino factory as well as 130 Km baseline for superbeam in presence of standard and 
non-standard interactions. For the oscillation channel considered by us the detection of muon is required. 
However, this is supposed to be easier than the detection of electron which is required in superbeam. This 
is an advantage in considering — )• signal in the neutrino factory. There are some studies on the 
performance of low enegy neutrino factories in the context of standard [5] and non-standard interactions [6] 
mainly for small 6*13. We have discussed the prospects of different baselines of neutrino factory considering 
large 613 as obtained from Daya Bay experiment. 

The paper is organized as follows: In section II we discuss Ve — >■ ly'^ oscillation probability and how the 
6 dependent and independent part varies with the variation of matter density for baselines L > O (2000) 
Km for standard and non-standard interactions. In section III we discuss the experimental set-ups and 
the assumptions in doing the numerical simulations using GLoBES. In section IV for standard interactions 
we have shown the discovery reach of CP violation and the precision measurement of CP violating phase 
6 in different long baselines and also for 130 Km baseline in superbeam. In presence of non-standard 
interactions we show the discovery reach of CP violation depending on various NSI true values and also 
discuss the discovery reach of NSIs depending on 6 true values. In section V we conclude with remarks on 
which baselines could be better for discovering CP violation after considering neutrino factory as well as 
superbeam for neutrino source. 



II. i/e Vu, OCILLATION PROBABILITIES WITH NSI 



Apart from Standard Model (SM) Lagrangian density we consider the following non-standard fermion- 
neutrino interaction in matter defined by the Lagrangian: 

C'^si = -2V2GFei^^[h^Pf][PisrLiyc.] (1) 

where P € {L,R), L = , R = ^^^^ ^ , / = e,u,d and are termed as non-standard interactions 

(NSIs) parameters signifying the deviation from SM interactions. Model dependent and independent bound 
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[7-9] are obtained for these NSI parameters. These NSI parameters can be reduced to the effective parameters 
and can be written as: 

where Uf is the fcrmion number density and is the electron number density. As these NSIs modifies the 
interactions with matter from the Standard Model interactions the effective mass matrix for the neutrinos 
are changed and as such there will be change in the oscillation probability of different flavor of neutrinos. 
Although NSIs could be present at the source of neutrinos, during the propagation of neutrinos and also 
during detection of neutrinos [10] but as those effects are expected to be smaller at the source and detector 
due to their stringent constraints [9], we consider the NSI effect during the propagation of neutrinos only. 
In section IV in numerical simulations we shall consider the model independent allowed range of real values 
of different NSIs as mentioned in reference [9] for earth like matter. 

In vacuum, flavor eigenstates may be related to mass eigenstates of neutrinos f, as 

\va>=^^Uai\vi> \ i = l,2, 3, (3) 

i 

where U is PMNS matrix [1] which depends on three mixing angles ^12, ^23 and ^13 and one CP violating 

phase (5. The Hamiltonian due to standard (Hsm) and non-standard interactions (-ffjvs/) of neutrinos 
interacting with matter during propagation can be written in the flavor basis as: 



H = Hsm + H^si 



(4) 



where 



Hsm 



2E 




A 





















0/ 



(5) 



i?NSi = A I e*^ e^^ e^r \ (6) 
In equations (5) and (6) 

where m, is the mass of the i-th neutrino and A is considered due to the interaction of neutrinos with matter 
in SM, G_F is the Fermi constant and is the electron number density of matter . £ee) , £erj Sjujui ^jur 
and Err are considered due to the non-standard interaction (NSIs) of neutrinos with matter. In equation 
(6), ( * ) denotes complex conjugation. In our numerical analysis we have considered the NSIs - Ee/i, ^er and 
EfiT to be real. 

To discuss the variation of 5 dependent and independent part in the oscillation probability and for that 
which baseline could be suitable for discovery reach of CP violation in presence of standard and non-standrad 
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interactions we present below the oscillation probability P^^^,^^ for long baseline (L > O 2000 Km). To get 
these expressions of probability we have followed the perturbation method adopted in references [10-13]. We 
shall present the oscillation probability upto order considering small NSI of the order of a. Considering 
the the matter effect parameter A in the leading order of perturbation and NSI parameters Sa0 of the order 
of a one obtains 
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where 
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+ 2| 
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For CP violation there is difTerence of probability in the neutrino oscillation and probability of antineutrino 
oscillation. One can relate the oscillation probabilities for antineutrinos to those probabilities given for 
neutrinos above by the following relation: 



Pa0 = Pa0{Scp -Sep, A -A). 



(10) 



In addition, we also have to replace e^^ with their complex conjugates, in order to deduce the oscillation 
probability for the antineutrino, if one considers non-standard interaction during propagation. 

To estimate the order of magnitude of 8 dependent and b independent part in the above two oscillation 
probability following Daya Bay result we shall consider sin^is ~ yja. For only SM interactions, (i.e e^^ — > 0) 
in above expressions of oscillation probabilities one finds that the b dependence occurs at order of o?!"^. 
This order of dependence with b remains same in the difference of neutrino oscillation probabilities and 
antineutrino oscillation probabilities (represented by AP later). However, the b independent part in AP 
(which could mimic CP violation) is at order a. This happens due to matter effect through A for SM as can 
be seen from above expressions. Apart from this difference in the order dependence by a^/^, the parameter A 
plays a very crucial role in determining the variation of b dependent (which varies as 1 j A) and b independent 
part (which varies as 1/(1 — A)). We have ignored the variation of trigonometric functions (which depend on 
A) as the variation is much slower for A due to the baselines in the range of 2000 to 5000 Km. For relatively 
shorter baselines the A value is smaller and the b dependent part becomes more pronounced. Apparently, it 
seems 2000 or 3000 Km baselines will be better than 4000 or 5000 Km baselines. However, it is important 
to note that for relatively shorter baselines the first oscillation peak will also correspond to relatively lower 
neutrino energy where the neutrino flux may not be that good. It is precisely for this reason that 2000 
Km (as found in our numerical simulation) is not that good so far the discovery reach of CP violation is 
concerned. 

However, when NSIs are also taken into account one can see that b dependence in AP could occur at the 
order of o?!'^ through 02 and 03 containing terms in (8) for NSIs of the order of a. We have checked that 
for slightly higher NSIs of the order of ^/a using perturbation method the same b dependent terms in AP 
appears with and 03 in the oscillation probability for long baseline as given in (8) and this slightly higher 
NSI makes these terms at the order of a. which could compete with the b independent part (which could 
mimic CP violation) in AP for long baseline as that is also at the order of a. So presence of slightly higher 
NSIs of order ^Ja present in and 03 improves the discovery reach of CP violation for longer baseline. As 
ai and 03 contains NSIs like e^^, and e^r it is expected that in presence of these NSIs the long baseline could 
provide a better discovery reach for CP violation. 
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III. NUMERICAL SIMULATION 



In this work wc have considered a neutrino factory set-up with four kinds of basehnes which are : 2000, 
3000, 4000 and 5000 km with 5.5, 7.1, 8.5 and 9 GeV parent muons respectively with 5 x 10^^ number of 
stored muons and anti-muons decays per year. The energy of parent muons have been chosen in such a way 
that the energy corresponding to the peak of flux matches with the neutrino energy correspond to the 
first oscillation peak of the probability of oscillation (fixing the probability for Standard Model 

interactions only) for different baselines. We have taken a magnetized totally active scintillator detector 
of mass 25 kt with threshold energy 1 GeV. In doing the analysis we consider a signal efficiency of 94% 
in appearance and 0.1% in disappearance channels. As signal we have taken the Uf. — )• oscilation and 
I'e — >■ oscillations. As backgrounds wc consider the neutral current events and and i/^ disappearance 
events. The systematic uncertainties of 2.5% and energy calibration error of 0.01% has been considered for 
both signal and for the background channels. The Gaussian energy resolution is considered to be O.ly/E . 
The numerical simulation has been done by using GLoBES [14]. For the experimental setup of CERN to 
Frejus baseline of 130 Km and for 2300 Km for neutrinos and antineutrinos coming from superbeam we have 
considered the flux and detector specifications as considered in reference [15]. For discovery reach of CP 
violation we have compared with 130 Km baseline because there the discovery reach is found to be better in 
comparison to other baselines. For NSI discovery we have compared with 2300 Km baseline as the discovery 
limit is also found to be better in comparison to other baselines for superbeam. 

We consider the true values of the neutrino oscillation parameters as lAml^l = 2.45 x 10~^ eV^, Am2i = 
7.64 X 10~^ cV^, 6i3 = 9°, 9i2 = 34.2° and 623 = 45°. Also in calculating the priors we consider an error of 
3% on 612, 0.005 on sin^20i3, 8% on 623, 4% on lAmg^l and 2.5% on Am^i- Also we consider an error of 
2% on matter density . 
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FIG. 1: Discovery reach of CP violation for only SM interactions for different baselines for normal (left panel) and 
inverted hierarchy (right panel). 



In the figure 1 we have shown the CP violation discovery reach for 2000 Km, 3000 Km, 4000 Km and 5000 
Km baselines for neutrino factory and 130 Km baseline for superbeam. It is found that 3000 Km baseline 



7 



for neutrino factory gives the better CP violation discovery reach than ah other basehnes considered for 
both the hierarchies of neutrino masses and for normal (inverted) hierarchy this discovery is possible over 
the 72 % (79%) of the allowed 6 values. However, 130 Km baselines for superbeam seems better than other 
baselines for neutrino factory as shown above. 




8 (true) 5 (true) 

FIG. 2: Precision of measurement of phase 5 for only SM interactions for different baselines for normal (left panel) 
and inverted hierarchy (right panel). 



In figure 2 we have shown at 5 cr confidence level the precision of measurement of phase S for different 
baselines as considered for CP violation discovery reach. From the figures one may find out the precision 
(Psitrue)) of sensitivity of measurement for any true value of 6 using the following expression for it: 

p 6{test){max) — 5{test)(min) 

S(true) 27J. 5{test){max) + 5{test){min) 

where S(test)(max) and 5{test){fnin) are the maximum and minimum 6(test) values respectively corre- 
sponding to certain true values. Here also it is found that for most of the 5(true) values the precision of 
measurement is found to be better for 3000 Km baselines particularly for normal hierarchy in comparison 
to other baselines. As for example for S (true) =0 the precision for normal(inverted) hierarchy is about 26 
% (18%) and for S (true) =tt/A the precision for normal(inverted) hierarchy is about 9 % (14%). Although 
for inverted hierarchy 4000 Km baseline is slightly better. 

In figure 3 we have shown at 5(j confidence level the discovery reach of CP violation in presence of NSIs 
Eee, Ee^i and for baseliucs 3000 Km and 5000 Km baselines for neutrino factory and 130 Km baseline for 
superbeam. 

In figure 4 we have shown at 5a confidence level the discovery reach of CP violation in presence of NSIs 
E/iM' ^tJ--r ^i^d for baselines 3000 Km and 5000 Km baselines for neutrino factory and 130 Km baseline for 
superbeam. In general, it is found that 3000 Km baseline in neutrino factory is somewhat better for smaller 
NSI values. However, for higher allowed NSI values 130 Km baseline for superbeam is found to be better 
in general. But for and Eer the CP violation discovery reach is much better for almost all allowed NSI 
values. The possible reason for this significant improvement is discussed in section II in the context of the 
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FIG. 3: Discovery limit of CP violation with NSI for both NH (left panel) and IH (right panel) at 5cr confidence 
levels for and 
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FIG. 4: Discovery limit of CP violation with NSI for both NH (left panel) and IH (right panel) at 5cr confidence 
levels for e^^, e^r and 
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oscillation probability. For 3000 Km baseline for e^e, Ce/j, Cer, Sfir and err the relatively better discovery 
reach of CP violation is possible at their NSI values -1, ±0.3, 0.5, -0.6, -0.3 and 0.5 respectively for both the 
hierarchies over 84%(82%), 97%(97%), 98%(98%), 69%(75%), 75%(75%) and 86%(83%) respectively of the 
allowed S values for normal (inverted) hierarchy . 

In figure 5 we have shown at 5a confidence level the discovery of NSIs e^e, ^en and e^r for baselines 3000 
Km and 5000 Km baselines for neutrino factory and 2300 Km baseline for superbeam. 

In figure 6 we have shown at 5a confidence level the discovery of NSIs e^^, e^r and Crr for baselines 3000 
Km and 5000 Km baselines for neutrino factory and 2300 Km baseline for superbeam. In general the NSI 
discovery reach is better for 3000 Km baseline with neutrino factory for different NSIs except (for which 
5000 Km baseline in neutrno factory is better) and e mur (for which 2300 Km baseline in superbeam is 
better). In case of 3000 Km baseline for See the NSI discovery is possible for its value above about ±0.1 
except S values around 225°(225°) and 80°(50°) for normal (inverted) hierarchy where the respective NSI 
values are slightly higher. For 3000 Km bcLseline for 6^^^,, Cqt^Cij,^^ Ci^r and Crr the discovery is possible below 
-0.02 and above 0.02 (below -0.02 and above 0.02), below -0.04 and above 0.02 (below -0.02 and above 0.05), 
below -0.18 and above 0.18 (below -0.18 and above 0.2), below -0.015 and above 0.015 (below -0.015 and 
above 0.015), below -0.18 and above 0.18(below -0.18 and above 0.18) respectively for normal (inverted) 
hierarchy. 



IV. CONCLUSION 



We have made a comparative study of different baselines 2000, 3000, 4000 and 5000 Km baselines for 
neutrino factory and 130 Km and 2300 Km baselines for superbeam in studying discovery reach of CP 
violation, precision of measurement of 5 for standard interactions. We have also made comparative study of 
baselines 3000, 5000 Km for neutrino factory with 130 Km for superbeam for CP violation discovery reach 
and with 2300 Km baseline for superbeam for NSI discovery limits. We have considered 130 Km baseline with 
super beam as it is found to give bettor discovery reach in comparison to other baselines for superbeam [15]. 
Similarly we have considered 2300 Km baseline with superbeam as it is found to give better NSI discovery 
limits in comparison to other baselines for superbeam [15]. In comparing baselines for neutrino factory with 
baselines for superbeam we have considered some conservative setups as mentioned in GLoBES. However, if 
the technology is further developed for superbeam [16] our conclusion could differ. 

Considering standard model interactions only we find that for CP violation discovery reach 3000 Km 
baseline for neutrino factory is somewhat better than other baselines with neutrino factory or with superbeam 
as found in figure 1. However, for precision of measurement of different values of 5 (true) both 3000 and 
4000 Km baselines for neutrino factory seem to be good as seen in figure 2. In presence of non-standard 
interactions over the entire allowed range and e^T for most of the allowed region 3000 Km baseline for 
neutrino factory is found to be good for CP violation discovery reach. The possible reason for this feature 
has been discussed in section II from the expression of — > oscillation probability. One advantage of 
neutrino factory also lies in the detector for identification of muon instead of electron which is required in 
superbeam. We find that the NSI discovery reach is better for 3000 Km baseline with neutrino factory for 
different NSIs except e^^ and e^^. So it seems that if we consider CP violation discovery reach, precision of 
measurement of 6 as well NSI discoveries 3000 Km baseline in neutrino factory performs better in most of 
the cases for various allowed range of NSI parameters. 
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